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Large spatial density fluctuations in high-energy heavy-ion coUisions can come from 
many sources: initial transverse density fluctuations, non-central collisions, phase tran- 
sitions, surface tension, and fragmentations. The common presence of some of these 
sources in high-energy heavy-ion collisions suggests that large scale density fluctuations 
may often occur. The detection of large density fluctuations by single-event Hanbury- 
Brown-Twiss interferometry in heavy-ion collisions will provide useful information on 
density fluctuations and the dynamics of heavy-ion collisions. 

1. Introduction 

Large spatial density fluctuations in high-energy heavy-ion collisions can come from 
many sources. Some of these sources are expected to be quite common in high- 
energy heavy-ion collisions and consequently large scale density fluctuations may 
often occur. It is of great interest to examine the sources of large density fluctuations 
and to fi nd ways for their detection so as to obtain useful information on the collision 
dynamicJiMSEl. 

Among the many different ways to study a quark-gluon plasma, intensity in- 
terferometry (HBT interferometry) is one of the experimental tools to examine the 

7 

space-time density distribution of the produced matter . The usual application of 
the HBT interferometry uses pairs of identical particles from many events (multi- 
event analysis) of similar general characteristics to study the (averaged) space-time 
configuration. Because of the difference in initial conditions in different collision 
events, different collisions will lead to different density fluctuations for some degrees 
of freedom while retaining similar fluctuations for some other degrees of freedom. 
An average over many events in a standard HBT analysis will lead to a distribution 
with suppressed fluctuations in event-dependent degrees of freedom but will retain 
those fluctuations in event-independent degrees of freedom. The proposed single- 
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event HBT analysis, if it can be carried out with sufficient accuracy, allows the 
examination of the full characteristics of large scale density fluctuations without 
suppression. The observation of HBT interferometry in atomic systems also make 
it interesting to examine the interference pattern in systems with large macroscopic 
density fluctuation&Sj^. 

If the statistics of identical bosons turns out to be insufficient for single-event 
HBT analyses using presently available detectors, it will be of interest to carry out 
few-event HBT analyses to examine event-independent density fluctuations. The 
few-event analysis will however require a good theoretical understanding of the 
single-event HBT interferometry. 

With the example of a granular density distribution to indicate the type of 
correlations one can get in a large scale density fluctuation, we shall show that the 
correlation function of a highly fluctuating density distribution is likewise highly 
fluctuating, with maxima and minima at locations which depend on the relative 
coordinates of the density centers. These local maxima and minima arise from the 
constructive and destructive interferences of identical bosons originating from two 
different centers. Their interference patterns therefore provide useful information 
on highly fluctuation density distributions in high-energy heavy-ion collisions. 



2. Sources of Large Density Fluctuations 



In many simulations of the heavy-ion collisions on an cvcnt-by-event basis, the ini- 
tial transverse energy density is far from being uniform. I t exh ibits large transverse 
density fluctuations with a large peak-to-valley A number of trans- 

verse density 'lumps' are clearly visible in t he i nitial transverse d ensi ty distribution 

or Fig. 4 of Ref. shown in Fig. 



m 



of the produced matter in Fig. 21 of Ref. 
1(a). 

Because of Lorentz contraction, the lon- 
gitudinal length of the produced matter is 
much smaller than the transverse length. Ini- 
tially as depicted in Fig. 1(a), density fluc- 
tuations in the transverse pla ne ma nifest 
as a number of transverse lumpj^^mi. The 
subsequent longitudinal expansion proceeds 
much faster than the transverse expansion, as 
pointed out by Landau and Belenkii - . The 
expanding matter develops into tubes [Fig. 
1(6)] and the density of matter in the tube 
decreases as a function of the proper time. 
The sausage instability leads to the break- 
up of the tube of matter into many approxi- 
mately equal spherical droplet depicted in 
Fig. 1(c), as in the break-up of toroidal liquid dropJ^^I, 



(a) Initial 



(b) Shortly afterwards 




(c) Later afterwards 



iniii 



riiin 



Fig. 1. The scenario in a central 
high-energy heavy-ion collision. 
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Another initial transverse density distribution occurs in non-central collisions of 
two approximately equal size nuclei. In these collisions, while the spectator matter 
proceeds forward, the participants collide and interact. The pion-emitting source 
produced by the participants may retain part of the initial momenta in different spa- 
tial regions. If the produced matter constituents retain a large fraction of the initial 
momenta, they can be approximately ideal- 
ized as two pion emitting clusters as shown 
in Fig. 2a. On the other hand, if a significant 
fraction of the participant matter in the over- 
lapping region is stopped while the more pe- 
ripheral constituents moves in their original 
directions, the pion emitting source can be 
approximated as consisting of three clusters 
as in Fig. 26. Interference of identical pions 
emitted from either two or three clusters will 
exhibit features associated with these spatial 
configurations. 




Fig. 2. Possible configurations of 
the pion-emitting source after a 
non-central collision. 



Another source of density fluctuation 
arises from the occurrence of first-order 

phase trarLsition, a s was first pointed out by WitterP^ and examined by many 
workerjmMia. According to such a description, incipient bubbles of low-density 
matter (hadrons) are formed in some local regions in a nearly static quark-gluon 
plasma, as the local temperature decreases below the phase transition temperature. 
The coalescence of these bubbles leads to bubbles of larger radii. As the system cools 
further, the fraction of matter in the low-density bubbles will become greater than 
the fraction of high-density matter (QGP), and the high-density matter will turn 
into lumps in the low-density matter background. The density of such a configura- 
tion is highly fluctuating. 

Still another source of density fluctuation arises from the effects of surface ten- 
sion during the evolution of the strongly-interacting matter. Because the QGP and 
the hadron matter are strongly interacting dense media, a surface tension arises 
at a boundary due to the presence of a strong interaction in the dense phase on 
one side of the boundary and the absence (or weakening) of the strong interaction 
with no density (or diminishing density) on the other side. This imbalance of the 
forces acting on different parts at the boundary leads to the surface tension, and 
the detail profile of the boundary may depend on the nature and the order of phase 
transition. Due to the presence of this surface tension, there may be the instabil- 
ity against surface shape changes that favors the formation of surfaces of smallest 
areas in different local regions. Such a tendency lead to large spatial density fluctu- 
ations. The end product may be the multi-fragmentation of the strongly-interacting 
matter. Mu lti-fragmentation has been observed in intermediate energy heavy-ion 
coUisionJI^. 
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3. Detection of Large Density Fluctuations by Single-Event HBT 
Interferometry 

In order to understand how the features of the correlation function may be affected 
by the presence of large density fluctuations, it is worth studying the correlation 
function of identical particles of momentum ki and ^2 with q = ki — k2 and P = 
{ki + ^2)72 for an idealized chaotic droplet sourcei^. There is a relation between 
the two-boson correlation function and the Fourier transform of the source density 
p{x] fci, ^2), which we assume to be independent of g, 



C7(q,P) = l 



(1) 



For brevity of notation, we shall leave the label P implicit and understood. 
We consider a density distribution of N droplets of the type 



N 

p{x) = A^pj{x), 



(2) 



where pj is the density distribution of the j/-th droplet, and A is a normalization 
constant such that the total density p[x) is normalized to unity as J dx p{x) — 1. 
We can consider the j-th droplet to have particle emission time centering at Tj , to 
be localized initially at Rj with standard deviations cij and Tj, and to move with 
a velocity Vj , 

Pj{t,v)^ ^ ^ j= . (3) 



27ra,)3 



2'KT, 



The correlation function is then 



C{q) 



N 

A^exp<; iq- 

.7=1 



(4) 



In order to get a clear insight into the most important features of the correlation 
function, wc consider the simple case where the density distributions of all droplets 
are the same so that aj — a and Tj=T for all j. In this simple case, A = 1/iV and 
the correlation function can be easily evaluated to give 



C{q) = 1 



7V2 



N 



(go-q-Vj) 



N 

E 

j,k=i 

j>k 



cos Aji, 



(6) 



Thus, the correlation function C{q) has maxima at Ajk ~ 2m:, and minima at 
Ajk ~ (2n — 1)tt, with n = 1, 2, 3, .... Numerical examples with = 4 and = 8 
static droplets are shown in Figs. 1 and 2 of Ref. {^\. The number of correlation 
function maxima for 8 droplets is greater than the number of maxima for 4 droplets. 
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The smaller the number of droplets, the greater will be the magnitude of fluctua- 
tions, as can be easily deduced from Eq. ([5]). Previously, we examine the two-pion 
interferometry for a granular source of quark-gluon plasma droplets and found that 
the granular model of the emitting source may provide an explanation to the RHIC 
HBT i?out / -Rside puzzle . We also investigate two-pion Bose- Einstein correlations 
of many droplets in single-event measurements. We find that the distribution of the 
fluctuation between correlation functions of the single- and mixed-events provide 
useful signals to detect the granular structure of the soured. 




(a) q, = 0, (b) q, = 0.05, (c) q, = 0.10 
Fig. 3. The correlation function for a two static droplet source. 

The case of iV = 2 presents an interesting example of the fluctuation of the 
HBT correlation functions. For the case of two static droplets of cr = 2 fm at an 
approximately touching configuration with a separation of 2\/5cr = 8.84 fm in the 
X direction. Figs. (3a) give the correlation functions C(q) at go = as a function 
for = 0, (36) for = 0.05 GeV/c, and (3c) for q^ = 0.10 GeV/c. All momenta 
in Fig. 3 are in units of GeV/c. In addition to the maximum at q = 0, there are 
prominent maxima and minima at q^ = 2tt/{Ri — -R2). For this case, the correlation 
function has minimum at q^ ~ 0.70 GeV/c and a maximum at q^ ^ 0.14 GeV/c. 
The correlation function decreases in magnitude as qz increases. 

In the case of two equal size droplets whose centers are moving with equal 
and opposite velocities vi = — V2 = ve^ along the beam direction (z-axis), the 
correlation function C'{q) is 

/2i2 2\2 22 22, 

C{q) = 1 + 2 \ C0sh{2T^ qaq^v) 

+ cos[q-i-R,-Ii2)-qo{Ti-T2)+qMTi-T2)]\. (7) 
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In this case, the correlation function C{q) is characterized by an asymmetric dis- 
tribution that is narrower in the Qz direction than in the transverse direction, in 
addition to a cosine-type modulation along the direction of the initial separation of 
the two droplets. The degree of the narrowing of the Qz distribution width depends 
on the of the emission time width r. We show in Fig. 4 the correlation function 
C{q) at go = for two droplets with the same geometrical dimensions as those in 
Fig. 3 but with Vz ^ 1 and t ~ a. The general features of the correlation function 
fluctuations are similar to the static case except that the correlation function falls 
faster as a function of Qz, as one can observe by comparing Figs. (3c) and (4c). 




Fig. 4. The correlation function for a two droplet source with Vz/c — 1 
and T = a. 

In a single-event HBT measurement, one can locate the reaction plane and 
designate the transverse axis in the reaction plane to be the x-axis. This a;-axis 
can be considered as the axis joining the centers of the two clusters initially. Then 
the correlation function will show prominent fluctuations in qx, in addition to the 
Gaussian distribution in Qz and qx, according to Eq. ([7]). 

If the statistics of single-event is insufHcient, one can perform an "aligned" 
multi-event HBT analysis by referring to the transverse axis in the reaction plane 
of each and every event as the x-axis. Then, fluctuations in qx in different events are 
event-independent and will be retained. A multi-event analysis with the properly 
aligned transverse axis will exhibit large fluctuations as those shown in Figs. 3 and 
4 and indicated by Eq. ([7]). 

On the other hand, one can use an "unaligned" multi-event HBT analysis in 
which we do not align the transverse axis on the reaction plane as the x-axis, 
Ri — R2 from different events can be considered randomly oriented in the transverse 
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plane. If the emission time Ti are the same, then the cosine term in Eq. ([7]) wiU be 
averaged to zero, and the correlation function becomes 

/2i 2 2\ 2 2 2 2 2 

C{q) = 1 + COsh{2T\oq,v). (8) 

One can extract the radii i?,;(IIBT) from this type of "non-aligned" multi-event 
HBT analysis at qo ~ 0, corresponding approximately to the inverse of the coeffi- 
cient of qf in the exponential factor in Eq. ([H]), and one obtains 

Rx ^ Ry = i?sido = and i?iong ~ VO^+T^U^, (9) 

i?iong/i?sidc ~ V^^T^/a. (10) 

This ratio is always greater than unity. Thus, while the single-event and the aligned 
multi-event analyses provide information on the fluctuation of the correlation func- 
tion, the unaligned multi-event analysis will give the i?iong/-Rsidc ratio to provide 
useful information on the emission time and the relative velocity of the two clusters. 

4. Discussions and Conclusions 

Density fluctuation occurs in high-energy heavy-ion collisions in many ways. It 
can arise from initial transverse density fluctuations as a result of the underlying 
granular structure of the colliding constituent nucleons. It can come from the non- 
central character of the collisions, leading to two or three clusters of pion emitting 
clusters. It can arise from the occurrence of first-order phase transitions. It can come 
from the surface tension effects which tend to form clusters of strongly-interacting 
matter during its evolution. The surface tension effects also lead to the multi- 
fragmentation of the produced matter into small clusters. Some of these density 
fiuctuation sources are expected to be commonly present and may be detected by 
single-event HBT interferometry. 

In a nearly head-on collision at very high energies, the number of identical pions 
is of the order of a few thousand. The number of observed identical pions is a 
small fraction of this number. For example, the number of identical pions detected 
in the STAR Collaboration in the most central Au-Au collisions at RHIC is of 
the order of a few hundred^. Although the number of pairs of identical pions in 
the event varies as n-^{nTr — l)/2, only a small fraction of these pairs have relative 
momenta small enough to be useful in a HBT analysis. Clearly, whether or not 
single-event HBT measurements can be carried out remains to be tested. 

If the statistics of identical bosons turns out to be insufficient for single-event 
analyses using the present detectors, there remains the future prospect of perform- 
ing single-event HBT analyses at RHIC with detector upgrades or with heavy-ion 
collisions at LHC. It will also be of great interest to carry out few-event HBT anal- 
yses by using pairs of pions in events with similar global characteristics. Because 
of the difference in initial conditions in different collision events, different collisions 
will lead to different density fiuctuations for some degrees of freedom while retain- 
ing similar fiuctuations for some other degrees of freedom. The use of few-event 
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HBT analyses may still retain those fluctuations from event-independent degrees of 
freedom. This may be particularly useful in the examination of two or three cluster 
configurations in non-central collisions by collecting events of similar impact pa- 
rameters and looking at pion pairs from these events. While the single-event and 
the aligned multi-event analyses provide information on the fluctuation of the cor- 
relation function, the multi-event analysis will give the i?iong/^side ratio which will 
provide useful information on the emission time and the relative velocity of the two 
clusters. 
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